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On 



The notion of vacuum fluctuations of the gravitational field plays important role in cosmology. The strong variable gravitational 
field of the very early Universe amplifies these fluctuations and transforms them into macroscopical cosmological perturbations. 
Since the underlying process is the parametric ampliflcation, the perturbations are placed in squeezed vacuum quantum states. 
It is possible that we already see the manifestations of these processes in certain cosmological observations. 
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I. ELECTRODYNAMICS AND 
GRAVIDYNAMICS 

The concept of the zero-point (vacuum) quantum fluc- 
tuations is best known in the context of theoretical and 
experimental studies of the electromagnetic field. Even 
in this case, there still remain delicate questions about 
the meaning of vacuum fluctuations and whether or not 
they produce, like other forms of matter, their own grav- 
itational field. Needless to say, vacuum fluctuations of 
the gravitational field itself seem to be even more enig- 
matic. The very quantization of the gravitational field is 
often considered as something mysterious. One will often 
be reminded that a fully consistent theory of quantum 
gravity is still a fundamental unsolved problem. Never- 
theless, there are some reasons to believe that we may 
be facing the manifestations of the gravitational vacuum 
fluctuations already now, in cosmological observations. 
To make our point clearer, we will start from the pre- 
sentation of the gravitational theory — the Einstein's 
general relativity — in the form similar to the theory of 
the electromagnetic field. 

The laws of electrostatics are expressed by the Poisson 
equation for the electrostatic potential 4>: 



Acj) = -Airp 



(1) 



where p is the charge density and A is the Laplace op- 
erator, A — + + On the route to electrody- 
namics, one replaces 1 electrostatic potential (?i by 4 elec- 
trodynamic potentials A^^ and the 1-component source p 
by the 4-component source j'^ (Greek indices p. u . . . run 
from to 3). One achieves the relativistic generalization 
of the theory by replacing the Laplace operator A by the 
D'Alambert operator □, □ = — + A, thus making 
the electromagnetic waves possible. As a result, one ar- 
rives at the equations of electrodynamics, the Maxwell 
equations: 



47r 

c 



or, in other notations, 



47r 
c 



(2) 



The Maxwell equations written in the form (2) assume 
that the gauge condition = is chosen. Otherwise, 
without fixing a gauge, the general form of the equations 
is 



47r 

^ ,Q ^ ,a — J 



(3) 



This form of equations makes the conservation law 
^ — automatically satisfied as a consequence of the 
equations. 

The 4- vector potential A^^ and the 4-vector current 
are defined in the Minkowski space-time 



dfj^ ^ ^rj.a ^ ^2 ^^2 _ ^^2 _ ^y2 



(4) 



they are functions of the Lorentzian coordinates 
t, X, y, z. If for some reason one needs to introduce 
arbitrary curvilinear coordinates in the Minkowski space- 
time, 



dcP' = dx^ dx^ 



(5) 



one is free to do this. Under arbitrary transformation of 
coordinates, the objects A^^ and transform as compo- 
nents of the 4- vectors. In Eq. (3), one will need to replace 
ordinary derivatives by covariant derivatives defined with 
respect to the flat space-time metric (5) written in curvi- 
linear coordinates. 

Being guided by the transition from electrostatics to 
electrodynamics, we can now try to guess how might look 
like the equations of gravidynamics, that is the relativis- 
tic field-theoretical generalization of the Newton law. 
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The gravistatics is described by the equation similar 
to Eq. (1): 

= 47rGp 

where (f> is the gravitational potential, and p is the mass 
density. Since the mass density p is only one component 
of the 10-component quantity, the energy-momentum 
tensor t^^ , one expects that the source term of the 
relativistic gravitational field will be the full energy- 
momentum tensor r^"^. Accordingly, one needs to re- 
place 1 Newtonian potential (/> by 10 gravidynamic po- 
tentials hi^" . The 10 components of h'^" comprise, simi- 
larly to T^'^, a symmetric second rank tensor. To make 
the theory relativistic and to allow for a finite velocity 
of propagation and waves, one needs to replace the el- 
liptic differential (Laplace) operator A by the hyperbolic 
differential (D'Alambert) operator □. And finally, be- 
ing aware of the universal character of the gravitational 
field, one expects that the energy-momentum tensor 
of the gravitational field itself should participate on the 
equal footing with the t^^" as the additional source term 
in the gravitational equations, thus making the equations 
nonlinear. One can guess that the equations of gravidy- 
namics should have the form 



ah"" = - 



167rG 



By analogy with the transition from electrostatics to elec- 
trodynamics, one may expect that this form of equations 
assumes that the gauge condition ft-^'^.iy = is already 
chosen. In arbitrary gauge, the general form of equations 
is likely to be 



IGttC / ,,,, ,,T,\ 



(6) 



This form of equations guarantees that the conservation 
law {t"" +t"'')^v = is a direct consequence of the equa- 
tions. It is obtained by taking derivatives from both sides 
of Eq. (6). 

The quantities ft-'^'^, r^'^, and t"" are assumed to be 
given in the Minkowski space-time (4). They are func- 
tions of the Lorentzian coordinates t, x, y, z. However, 
one is free to introduce arbitrary curvilinear coordinates 
in the Minkowski space-time and to write the flat space- 
time metric in the form (5). Under arbitrary coordinate 
transformation, the quantities /i'"', t"'^ , t"" transform as 
symmetric second rank tensors. To write Eq. (6) in arbi- 
trary curvilinear coordinates, one needs to replace ordi- 
nary derivatives (coma's) by covariant derivatives (semi- 
colon's) defined with respect to the metric (5). As a 
result, one obtains the equations 



167rG 



Using the variational principle, Eqs. (7) can be derived 
from the action S consisting of the gravitational part 
and the matter part (including interaction with gravity) 



5 = 53+5™, 5^ 



167rG 



(7) 



Specific Lagrangians L^, determine, of course, the 
entire content of the theory, but we will not need 
their explicit form in our brief discussion. The energy- 
momentum tensors f^" , t"" are defined in a traditional 
field-theoretical manner, as variational derivatives of La- 
grangians with respect to the metric tensor 7^^. The 
energy-momentum tensor of the gravitational field is 
tf_i,i, = — the energy- momentum tensor of 

the matter soiirces (including their interaction with grav- 
ity) is r^j,y = Like in the case of electrodynam- 
ics, gravidynamics admits (in addition to the freedom 
of making arbitrary transformations of coordinates) the 
gauge freedom, that is a possibility of replacing the par- 
ticipating dynamical variables by new ones constructed 
from the old ones according to a certain prescription 
(gauge transformations). One can use this freedom for 
imposing gauge conditions, for instance, h""-^,^ = 0. 

The point of the discussion above is that Eqs. (7), (or 
Eqs. (6) if we use the simplest, Lorentzian coordinates) 
which we envisaged on physical grounds, are exactly the 
Einstein equations, although written in a somewhat un- 
familiar fashion. We have outlined the field-theoretical 
formulation of the general relativity, in contrast to the 
usual geometrical formulation. The link between the two 
formulations is realized by introducing the new functions 
^—g g"" = -x/— 7(7^^ + h"") and interpreting the func- 
tions gijiu as the metric tensor of the curved space-time: 
ds^ = gfivdx'^ dx"^. The matter energy-momentum tensor 
T'"^ participating in the geometrical Einstein equations 
is defined as the variational derivative of the with 
respect to what is now metric, that is g"'^ . By using 
(7''^ and T^'', one can transform Eqs. (7) to the familiar 
geometrical representation of the Einstein equations: 

The economy and beauty of the geometrical formulation, 
in which a single object g^,^ plays the double role of 
the metric tensor and the functions describing gravity, 
is accompanied by the loss of the gravitational energy- 
momentum tensor t^'^ in exchange to ugly constructions 
called pseudotensors. 

The two formulations of the general relativity are es- 
sentially two sides of the same theory. The use of one 
or another framework becomes a matter of taste. For 
our purposes, it will be convenient to think of the gen- 
eral relativity in terms of the field-theoretical approach, 
as of gravidynamics, similar in many aspects to electro- 
dynamics. We will discuss various cosmological models, 
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but from the field-theoretical position, they are simply 
specific gravitational fields. For instance, the averaged 
structure of the observed part of the Universe is well 
approximated by a homogeneous isotropic cosmological 
model which can be described, in geometrical language, 
as curved space-time with the metric tensor 

ds'^ = g^^ dx^ dx"" = c^dt^ - a^{t){dx'^ + dy^ + dz^) . 

(8) 

The function a{t) is called cosmological scale factor. This 
function shows how the distance between the remote cos- 
mological objects varies with time. We can also think of 
this cosmological model as of a gravitational field deter- 
mined by the gravitational potentials 

hoo = a^{t) - 1 , /ill = h22 = /133 = 1 - a,{t) , 

(9) 

specified in flat space-time with the metric tensor (4). 
From this position, the function a{t) describes the 
strength and variability of the gravitational field. In 
our further discussion this function will play the role of 
the pump field acting upon quantized perturbations of 
the cosmological gravitational field and whatever matter 
sources present in our model. 



II. AMPLIFICATION OF CLASSICAL 
COSMOLOGICAL PERTURBATIONS 

We happened to live in approximately homogeneous 
and isotropic Universe. The observed volume of space 
includes immense number of galaxies and various sorts of 
radiation. On scales of galaxies and their clusters, matter 
is concentrated in well discernible systems. However, on 
larger scales accessible for observations, the distributions 
of matter, its velocity, and accompanying gravitational 
field do not show any significantly preferred places or 
directions. The photons of the microwave background 
radiation traveled to us for 10 billion years and, yet, the 
temperature of the radiation visible in different direc- 
tions on the sky is remarkably the same. A priori, one 
would not be very surprised if one side of the observed 
sky contained, say, 20% more matter than the other. For 
instance, we know that our own Galaxy is a flattened 
system of stars in which the Sun is located not very far 
from the edge of the disk. And there would be noth- 
ing wrong in having temperature of the microwave back- 
ground radiation in one direction on the sky, say, 20% 
higher or lower than in others. However, this is not the 
case. On large scales, the deviations from homogeneity 
and isotropy are small. Most convincingly, this is demon- 
strated by the fact that the measured large angular scale 
anisotropics of the microwave background radiation tem- 
perature ST/T have the level of only 10"^ - lO"''. So, 
when we work with homogeneous isotropic cosmological 



models plus small perturbations, this is not just a mathe- 
matical simplification, this is a reflection of observational 
data about the real world. 

There are three ways in which one can perturb a homo- 
geneous and isotropic distribution of matter and fields. 
First, one can compress matter in various places per- 
turbing its mass density, velocity and the accompanying 
gravitational field. This is called cosmological density 
perturbations. Second, one can provide the elements of 
matter with small rotational velocities (without perturb- 
ing the mass density) which will also be accompanied by 
perturbations of gravitational field. This is called cos- 
mological rotational perturbations. And third, one can 
perturb gravitational field itself. This is called cosmolog- 
ical gravitational waves. As we know, the equations of 
gravidynamics allow gravitational waves. In contrast to 
gravitational waves emitted by localized sources, we will 
be dealing with waves as excitations of gravitational field 
in the entire Universe. 

Let us concentrate on gravitational waves — the entity 
analogous to electromagnetic waves. 

The gravitational field of an idealized, strictly homoge- 
neous and isotropic universe, is described by the poten- 
tials (9). We want to consider a more realistic universe, 
which includes gravity-wave perturbations. The poten- 
tials can now be written in the form ft,^'' = +h'"''^^^ 
where h'"''^°^ is the zero-order approximation (9) and 
JiiJ-^W are the gravity-wave corrections. Since we expect 
the perturbations to be small, we can neglect quadratic 
and higher-order terms of h'^"^^^ in Eq. (6). By using the 
available gauge freedom we can eliminate the components 
^oo(i)^ ^oi(i) (i^g^^i^ indices i,j, ... run from 1 to 3). The 
remaining components /i'^^^^ are functions of t, x, y, z 
and can be Fourier decomposed over spatial harmonics 
gjnx^ g-mx^ where n = {n^,n'^,n^) is arbitrary wavevec- 

tor and the wavenumber n is n = 

Let us write down one of the Fourier components: 



(10) 



The two polarization tensors Pij (s 
conditions 



1, 2) satisfy the 



0, Pi 



0, 



PijP 



2(5. 



Pij{-n)=Pij{n). 



Gravitational waves have two independent polarization 
states. 

We should now derive the equation governing the time 

s 

dependent complex amplitudes Mn(i)- We should use 
Eqs. (6) retaining in them only terms linear in hij'^^K 
The nonlinear right-hand side of the equations will pro- 
vide us with cross terms, products of the zero-order terms 
g^jjj ^Yie first-order terms hij'^^\ These cross terms 
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signify the interaction of the gravitational waves with the 
pump field. As a result, we obtain the equation 



fin + n fin 



-^-^ - + - + 



a . 

a 



(11) 



Both polarization amplitudes satisfy this equation. 

We can now introduce a new time parameter r] accord- 
ing to the relationship cAt = adrj and denote d/dr] by a 
prime. Equation (11) acquires the form 



a 
a 



Mn = 



(12) 



We have derived this equation using the field-theoretical 
approach, but this is exactly the same equation which 
could be derived (and was in fact originally derived) using 
the geometrical approach. In terms of ?7-time the metric 
(8) is manifestly conformally-Minkowskian: 

ds^ = 0^(7?) (rf??^ - dx^ - dy^ - dz^) . 

The form of the equation (12) explains why certain 
phenomena in cosmology and quantum optics are linked 
by common underlying physics. Eq. (12) can be inter- 
preted as equation for a parametrically excited oscillator. 
The term — a" /a can be considered as square of the 
time dependent frequency of the oscillator (even though 
in the most interesting regime this term is negative). This 
equation can also be treated as a kind of the Schrodinger 
equation for a particle with energy in the presence of 
the potential barrier U{ri) = a" /a. Here we will mostly 
discuss gravitational waves, but other cosmologic;al per- 
turbations (rotational perturbations and density pertur- 
bations) are described by equations of similar nature and 
many portions of our discussion apply to them as well. 

Let us show how the gravitational pump field can am- 
plify classical gravitational waves, leading essentially to 
the production of standing gravitational waves. We will 
illustrate this process on a specific example. This exam- 
ple is a currently popular cosmological model which may 
have some relevance to the real Universe. 

Let us assume that the function a{ri) has the following 
behavior 



-oo < T] <r]i 



(13) 



where _ff is a constant, the c/H has dimension of length. 
The a(ry) and its first time derivative a'(ry) are continuous 
functions at the point rj ~ rji. In cosmology, function (13) 
describes the so called De Sitter era with the scale factor 
ai{ri) followed, after the moment of time r] = rji, rji < 
0, by the radiation-dominated era with the scale factor 
a^{ri). The potential barrier U{rf) in Eq. (12) vanishes 
asymptotically for rj ^ — oo, reaches the maximum at 
7] = r]i, and drops to zero after the moment r] = rii. 



Equation (12) admits simple exact solutions in the 
both eras. For rj <r]i, one has 



Hi = A 



cos(n?7 + (f)) — sm{nr] + ^) 



(14) 



where A and </>, < </> < 27r are arbitrary constants. In 
the asymptotic regime r/ — > — oo, we have usual harmonic 
oscillations 



Hi Pi A cos(n?7 -|- (p) 
For T] >rii, one has 

fie = B sin{nr] + x) 



(15) 



(16) 



where B and x are also arbitrary constants. 

Solutions (14), (16) and their first time derivatives 
should join continuously at 77 = 771. This allows us to 
express the final amplitude B and phase x in terms of 
the initial amplitude A and phase (j>. It is convenient to 
introduce the notation x = niji . Then we have 



B\c^) 
where tan 2a 

ta.nx{<t>) = ^ 

where 
tan/3i 



A^ 

2x4 



[1 + 2x* - cos(20 + 2x- 2a)] 



(17) 



2x 
1-2X-' 



, and 



' 1 -h 2a;4 - cos(2a; + 2a) sin(^ -|- (3i ) 
l + 2x'^ + cos(2a; + 2a) sin(^T^ 



(18) 



sin x — X 
sin X cos X — X 



tan/32 



sin X cos X — X 



The parametric amplification is a phase sensitive pro- 
cess, the numerical value of B depends on the initial 
phase (j). However, the root mean squared (r.m.s.) fi- 
nal amplitude is always larger than A, that is a "typical" 
oscillation will always be amplified. To show this, inte- 
grate B^{(p) by (j) from to 27r, divide the result by 2n, 
and take the square root of this number. You will get 



Since \x\ <C 1 for the frequencies n of our interest, we 
obtain 

, A 

We can also calculate the variance of the amplitude and 
find out that the variance is large. Namely, 



54 - (^2)2 ^ ]_ (^2)2 ^ 

A^ ~ 2 ~A^ ~ 8^ 



» 1 
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that is the dispersion of the quadratic amphtudc is of 
the same order as the mean value of this quantity. 

Under the condition |a;| <§; 1, the final phase Xi as a 
function of the initial phase </>, will be concentrated near 
two values of x separated by tt. From Eq. (18) we have 
approximately 



h\f w MB\pij cos ( nr] + 



tanx 



-2x 



sin((^ + /3i) 
sm{(j) + P2) 



The function tanx stays constant and almost zero for all 
(j) except for the values of 4> near (j) = 0, it, 2tt, . . ., where 
tan X diverges and changes sign. Almost all initial phases 
4> from the interval between and tt produce ( "squeeze" 
to) the final phase Xi ~ 0) and almost all initial phases (f) 
from the interval between tt and 2-k produce ("squeeze" 
to) the final phase X2 ~ tt. 

So far, we have discussed the time dependence of the 
Fourier component (10). Let us now take into account 
the spatial dependence too. Consider, for simplicity, only 
one polarization state and write 



(19) 



At the early stage r] — > —00, take, in agreement with 

Eq. (15), 



fj, ~ A e 



(20) 



where A = \A\e The component (19) can be rear- 
ranged to read 

h^j^ « 2\A\pij cos(nJ7 - nx + 0) , 

that is we have initially a traveling (not standing) wave 
propagating in the direction n. In the limit rj +00, 
the solution (20) will have the form 



/X « e-'"" + B2 e 



inrj 



(21) 



Since we are discussing the asymptotics of the same so- 



lution, the Wronskian /i/x*' — 
same constant for 77 —> —00 and 77 

\Af = \B^f - IB^f 



const should be the 
> +00. This gives: 

(22) 



^From the point of view of physics, this condition ex- 
presses the conservation law of the linear momentum. As 
we know, a "typical" fin al am plitude is mu ch lar ger than 
the initial one, that is, jsTp > and \B^ > \A\'^ . 
This means, that, as a rule, the number on the left hand 
side of Eq. (22) is a result of almost complete cancella- 
tion of the two much bigger numbers on the right hand 
side of that equation. This allows us to neglect |Ap and 
to write 



« \B\e- 



-iipi 



B, 



IBle^'f'^ 



(23) 



Using these values of iJi, B2 in Eq. (21), we can rearrange 
Eq. (19) to read, in the limit r] — > -|-oo. 



COS nx -|- 



V'2 - -01 



(24) 

This formula demonstrates explicitly that, after the 
amplification process, we are dealing with a standing 
wave. Since we have omitted the term \A\^ in Eq. (22), 
neglecting in this way the initial momentum, expression 
(24) is a strict standing wave. Acting more accurately, 
we would get an almost standing wave. The amplified 
initial traveling wave combines with the generated wave 
traveling in the opposite direction to form an (almost) 
standing wave. If we have had a strict standing wave ini- 
tially, it would be amplified (except in the very narrow 
interval of inappropriate initial phases) maintaining its 
standing wave pattern. 

We have prepared ourselves for quantum-mechanical 
treatment of this problem. Classical process of ampli- 
fication requires a nonzero initial amplitude. Other- 
wise, the final amplitude will also be zero. Quantum- 
mechanically, we do not need initial amplitude, or, we 
can say, some minimal level of the amplitude is provided 
by the zero-point quantum fluctuations. Since we can 
think of a vacuum state as of a set of fluctuations with 
all phases, we can expect that the amplification process 
will necessarily take place. Quantum-mechanically, this 
process describes creation of strongly correlated pairs of 
waves (particles) by the external pump field. The already 
briefly discussed distributions of the final amplitudes and 
phases (assuming that the initial phase </> is distributed 
randomly and evenly in the interval — 27r) are related 
to the phenomenon of squeezing. 



III. QUANTIZED COSMOLOGICAL 
PERTURBATIONS AND SQUEEZING 

It is now clear what we are going to quantize — the 
fluctuating part /i^^(i) of the gravitational field. If you 
think of gravitation as of geometry, you may be lead to 
the idea of quantizing, emitting and absorbing, pieces of 
space and intervals of time. If you think of gravity as of 
a field, you are lead to the idea of quantizing the field in 
no more mysterious fashion (even if more difficult tech- 
nically) than quantization of the electromagnetic field. 
There is not much sense in quantizing the pump field 
h^^^^^\ we are not going to discuss such issues as creation 
of universes, we are going to discuss creation of excita- 
tions in the Universe. Our theory of small excitations is 
essentially linear. The fundamental nonlinearity of gravi- 
dynamics provides us with the parametric interaction of 
the quantized perturbations with the pump filed. In con- 
trast to quantum optics, we do not need any intervening 
nonlinear optical medium in order to couple the quan- 
tized part of the field with the pump field, gravity is 
inherently nonlinear. 

Let us write the gravity- wave field operator in the gen- 
eral form 



5 



C 1 /"OO ^ 



c„(77)e-- + c:(7?)e- 



(25) 



Being armed with the gravitational energy-momentum 
tensor t'^'^ , we can find the normahzation factor C. In 
the traditional manner, we require that the vacuum en- 
ergy of each mode be equal to a "half of the quantum" . 
From this requirement wo find C — VIStt I pi where Ipi 
is the Planck length, Ipi = (Gh/c^)"! . The gravitational 
constant G and velocity of light c participate in the def- 
inition of t^" . The Planck constant fi enters because of 
our normalization. All the three fundamental constants 
combine in Ipi and appear in G, explicitly demonstrating 
that we arc working with quantum gravity. 

The creation and annihilation operators Cn{r]), 0^,(77) 
for each polarization state are governed by the Heisen- 
berg equations of motion 



dr] 



= -i[cn{r]),H] 



dr] 



(26) 



where H is the Hamiltonian for the modes involved. This 
Hamiltonian follows from the general Lagrangian of our 
theory, truncated at the appropriate order. Combining 
those formulas with the understanding of what kind of 
processes we are dealing with, one can derive 

H = ncicn + ncl„c_„ + 2a{r])cl,c'_n + 2(7* (r?)cnC_n . 

The coupling function a{r]) is determined by the pump 
field: a{n) = 1^. 

The form of the Hamiltonian dictates the form of the 
solution (Bogoliubov transformation) to Eq. (26): 



C„(r?) = U„(??)Cn(0) + Vn{v)c^_n{0) , 

ci{v)=K{v)ci{0) + v:{v)c.^{0) . 



(27) 



The initial values Cn(0), cj,(0) of the operators are taken 

at the very early times, 77 — !■ — 00, long before the interac- 
tion became effective. The two complex functions Un{r]), 
Vnif]) satisfy the equations 



. dUn . a' ^ 

I—— = nUn + i—v„ 
drj a 

and the condition 



I—— = nvn + i—u^ (28) 
drj a 



1 



(29) 



which guarantees that the commutator relationship 
[cn(?7), cjj,(r?)] = (5^(n — n') is satisfied for all times. The 

condition (29) is automatically fulfilled if one introduces 
the (squeeze) parameters r, 0, 4> according to the defini- 
tion 



u = e cosh r 



V = e-^f^-^-?^) sinh r . 



In our problem, these parameters arc functions of time. 

Using Eqs. (27) one can present the field operator (25) 
in the form 



G f°° ^ s 

^ij ^) = (2^)3/2 J_ '^^n 



hnivV^ c„(0) + /i*„(r;)e-^"^ c„(0) 



where the functions /in('?) are 
1 



1 



2n 



(30) 



hniri) = 



a{r]) 



+ <*(??)] 



a(ry) 



Mn(??) 



^From equations (28) for Un{ri), ««(??) we recover classical 

s 

equations of motion (12) for the functions Mn(r?). 

We have reduced our problem to the problem stud- 
ied also in other areas of physics. By cosmological rea- 
sons, we make the simplest assumption about the initial 
state, namely that each mode of the field was initially in 

the vacuum state |0n) defined by Cn(0)|On) = 0. In the 
Schrodinger picture, the initial vacuum state |0„)|0_„) 
transforms into a two-mode squeezed vacuum quantum 

state. To see this formally, wc can note that the Bogoli- 
ubov transformation (27) can be presented as 

Cniv) = RS c„(0)S'^i?^ , ci{T]) = RS ci{0)S'<R'' 



where 

S{r, (j)) = exp r^e 



^c„(0)c_„(0)-e2^<^c],(0)cL„(0)) 



(31) 



is the two-mode squeeze operator, and 

R{e) = exp[-z0(c],(O)c„(O) + cl„(0)c_„(0)) 

is the two-mode rotation operator. In our case, the two 
modes are two waves traveling in the opposite directions. 
Creation of pairs of correlated waves (particles) corre- 
sponds to classical picture of generation of almost stand- 
ing waves. 

We can also work with the one-mode squeezed states 
and associated standing waves. To reach this goal, intro- 

s « 

duce two new sorts of operators, and hn, according to 

the rule 

s 1 fs s \ s 1 /s 

In terms of these operators, the field (25) can be written 



as 



^ -j^ poo 2 ^ ^ 

' V2a{rj){2n)i J-00 t[ ' V^^ 
a„(??) + a„(?7)j cosnx -|- [bniv) + bn iv)j sinnx 
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demonstrating the standing wave decomposition. Ac- 
cordingly, the two-mode squeeze operator (31) factorizes 
into a product of two identical one-mode squeeze opera- 
tors: 



^r(e-2'^ a„(0)a„(0) - e^'^ aj. (0)at (0)) 

\r{e-^'^ 6„(0)&„(0)-e2^<^ biml{Q) 



S{r, <p) = exp 
X exp 



The "particles" described by and bn operators do al- 
ready have zero linear momentum. 

We will work in the Heiscnberg picture, in which 
the initial vacuum state does not change. It immedi- 
ately follows from the representation (30) that the mean 
quantum-mechanical value of the field hij [rj, x) is zero 
in every spatial point and at every moment of time: 
(0| /lij (ry, x) 1 0) =0. As might have been expected, the 
variance of the field is not zero and does depend on time. 
Using a product of two expressions (30) one can show 
that 



(0|%(77,x)/i^^(,7,x)|0) = 



2 /-oo 

/ 



hn{Tl) 



dn . 



(32) 



The expression under the integral is usually called the 
power spectrum of the field (in this case, the power spec- 
trum of the quantity hijh^^): 



(J2 

P{n) = —^n 



(33) 



lu K, 10^* cm and beyond. Correspondingly, the squeeze 
parameter r ranges from r w 1 at the short-wavelength 
end up to r w 100 and more for extremely long waves. 
The r.m.s. numerical values of the field are sufficiently 
high to believe that the quantum-mechanically generated 
cosmological perturbations may provide explanation to 
some of cosmological observations. For instance, if the 
length parameter cjH in the simplest model (13) is 5 or 6 
orders of magnitude larger than /pj, the r.m.s. amplitude 
of the generated gravity- wave field is appropriate to pro- 
duce the large-angular-scale variations in the microwave 
background radiation at the level jT k. 10~^ - 10~^. 

The two most important features of strongly squeezed 
vacuum states — small variances of phase and large vari- 
ances of amplitude — are responsible for the character- 
istic features of cosmological perturbations themselves 
and of the microwave background anisotropics caused by 
them. The small variances of phase and the associated 
standing wave pattern of the field make the spectrum (33) 
an oscillating function of frequency n for a fixed moment 
of time ry. The possibility of seeing these oscillations di- 
rectly or indirectly is one of the observational challenges 
of cosmology. On the other hand, large variances of am- 
plitude should lead to large statistical deviations of the 
observable quantities from their expected values. One 
may recall that, in a squeezed vacuum state, the disper- 
sion of the number operator A'' is of the same order as 
the mean value of A^: (A'^) - {Nf = 2{N){{N) + 1). In 
the next section we will focus on the consequences of the 
large variances of amplitude for statistical properties of 
the microwave background anisotropies. 



The integral (32) is always divergent in the ultra-violet 
limit n oo, and sometimes divergent in the infra-red 
limit n ^ 0. The ultra-violet divergence is caused by 
the high-frequency waves (above the barrier in Eq. (12)) 
which were not affected by the amplification process 
and which survived from 77 ^ — oo to 77 ^ -|-oo. The 
renormalization of Eq. (32) (subtraction of the "half 
of the quantum" contributions from the final expres- 
sion) eliminates the divergence and this part of the spec- 
trum entirely, without changing substantially the lower- 
frequency part of the spectrum where significant ampli- 
fication took place. As for the infra-red divergence, it 
arises if one assumes that the amplification process had 
been taking place for infinitely long time, which does not 
happen in practice. So, in practice, the integral (32) has 
natural high-frequency and low-frequency cutoffs, and 
does not require renormalization. 

For cosmological applications it is important to know 
the r.m.s. values of the field, its statistical properties, 
shape of the spectrum. Of course, they all are determined 
by the strength and variability of the pump field, and 
by the very nature of squeezed vacuum quantum states. 
In realistic cosmological models, the generated spectrum 
extends from wavelengths of order of a few meters up 
to the wavelenghts comparable with the Hubble radius 



IV. STATISTICS OF THE MICROWAVE 
BACKGROUND ANISOTROPIES CAUSED BY 
SQUEEZED COSMOLOGICAL PERTURBATIONS 

Why do we think that the quantum gravity processes 
in the very early Universe may have something to do with 
current observations? 

The first point is that the mechanism studied above 
is universal, it does not require any hypothesis beyond 
the validity of general relativity and basic principles of 
quantum field theory. One can believe that, in one or an- 
other amount, the quantum-mechanically generated per- 
turbations should be present. We can not make an exact 
prediction about the perturbations, because the pump 
field (the scale factor of the very early Universe) is not 
known. (And if it were known, there would be not too 
much to learn by testing our predictions.) But even the 
absence of perturbations on a certain level tells us some- 
thing positive. By comparing theoretical predictions of 
various cosmological models with observations we restrict 
the possible behavior of the very early Universe rather 
than refine hypotheses on the generating mechanism. 

The second point is that the measured large-angular- 
scale anisotropies are produced by cosmological pertur- 
bations with extremely long wavelengths. It is difficult 
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to devise an alternative mechanism of their generation, 
especially if we do not want to make many additional as- 
sumptions. At the same time, the quantum-mechanical 
generation suggests itself as natural and unavoidable. 

If the observed large-angular-scale anisotropics are in- 
deed caused by cosmological perturbations of quantum- 
mechanical origin (further theoretical study shows that 
they should be predominantly gravitational waves), the 
anisotropics should carry signatures of the underlying 
squeezed vacuum quantum states. One possibility to 
prove (or disprove) the quantum-mechanical origin of the 
perturbations is to use their statistical properties. 

Let us first discuss how anisotropics arise in the pres- 
ence of the gravitational field perturbations, regardless 
of what is the origin of perturbations. In a strictly ho- 
mogeneous and isotropic universe, the temperature of 
the microwave background radiation seen in all direc- 
tions on the sky would be the same, T. The field (9) is 
independent of spatial coordinates and whatever its ef- 
fect on propagating photons is, it docs not depend on 
direction. In presence of additional gravity field poten- 
tials hij{r],x), the frequency of photons (and, hence, the 
temperature of the Planck distribution of the microwave 
radiation) changes along their paths. Even if the pho- 
tons were emitted with exactly the same frequency (de- 
termined by the physical conditions of local equilibrium 
in points of their release) they will arrive to us from dif- 
ferent directions with different frequencies, because they 
have traversed different paths in the field hij{r],x). In a 
given direction on the sky, characterized by a unit vector 
e = (e^,e^,e"^), the actual measured temperature differs 
from its unperturbed value T. Using either geometrical 
or field-theoretical approach, one can show that the vari- 
ation of temperature in a given direction e is expressed 
by the formula 

ST 1 /"""i dh ■ ■ ■ 

where dhij/drj is taken along the integration path = 
e'w, rj = rjj^ — w, from the event of reception w = to 
the event of emission w = wi = rjn — ije- 

In our case, the field hij entering Eq. (34) is a quantum- 
mechanical operator (30). Hence, the ST/T does also 
become an operator. Having defined in this way an ob- 
servable and knowing the quantum state of the field |0), 
we can compute various expectation values. 

Obviously, the mean value of ST/T in every direction 
on the sky is zero: 

(0|f^(e)|0) = . 

The product of two operators ^(ei) and ^(62) for two 
different directions ei and 62 is a new operator. The 
mean value of this operator is called the angular correla- 
tion function: 

i^(ei,e2) = (0|^(ei)^(e2)|0) . (35) 



One can calculate this quantity by manipulating with the 
product of two expressions (34) and performing integra- 
tion over angular variables in the 3-dimonsional wavevec- 
tor n space. One can show, first, that the quantity 
K{ei,e2) depends only on the separation angle 5 between 
the vectors ei, 62, not vectors themselves, and, second, 
that the quantity i^T (61,62) takes on the form 

00 

K{eue2) = K{6) = l%i ^ KiPi{cosS) . (36) 

^— ^min 

The coefficient Zp; is taken from C^, other numerical 
coefficients are included in if;. The P;(cos(5) are the 
Legendre polynomials. The lowest multipole imin follows 
from the theory and, in case of gravitational waves, it is 
^min = 2. The quantities Ki are calculable from a given 
cosmological model plus perturbations, and they even- 
tually depend on the pump field. For a zero separation 
angle (5 = 0, one obtains the variance of ^(e): 

{0\ — {e) — {em=K{0) = l%, J2 • 

i = imin 

The numerical level of ST/T is basically determined by 
the r.m.s. numerical level of the generated amplitudes of 
the perturbations. Large variances of the amplitude will 
be reflected in large statistical deviations of the measured 
values of the observable ^(61)^(62) from its expecta- 
tion value (36). Let us calculate the variance V^(ei, 62) of 
the operator v = ^(61)^(62). By definition, we have 

y(6i,62) = (0|v'|0)-[(0|v|0)]' . 

This calculation requires us to deal with the product of 
4 expressions (34). However, this calculation is easy to 
perform. One can show that, first, the quantity V(6i, 62) 
depends only on the separation angle S, and, second, its 
concrete value is 

V(ei,e2) = V{S) = K^(S) + K^O) . 

As could be expected, the standard deviation 

a{S) = ^/V{&j = ^/K\S)+K\0) 

is very big. Even for those separation angles for which 
the K(S) vanishes, the standard deviation of the vari- 
able V is as big as the mean value of this variable. If 
the anisotropics are indeed caused by the squeezed cos- 
mological perturbations, this statistical feature must be 
present in the observational data. The problem is how 
to recognize it having only a single map of the sky. The 
measured anisotropics are essentially a single function 
of angular coordinates which assigns different numerical 
values of the temperature to different directions on the 
sky. 

A complication, specific for cosmology, is that we will 
always be dealing with only one realization of our quan- 
tum state, or, in other words, with only one experiment. 
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only one map. The calculated expectation values im- 
ply averaging over many realizations, many experiments, 
"many universes", but we have access to only one of 
them. The meaningful comparison of the theory and ex- 
periment will require a kind of ergodic hypothesis allow- 
ing us to replace the quantum-mechanical expectation 
value of, say, observable v by the procedure of averaging 
the measured ^(ei)^(e2) over many pairs of directions 
ei, 82 on a single map. 

Even if it is known in advance that the anisotropics 
are induced by squeezed perturbations, the extraction 
of cosmological information from the observations is not 
easy. The large statistical uncertainties discussed above 
make the discrimination of competing cosmological mod- 
els (competing pump fields) a much more difficult prob- 
lem than we wanted it to be. It is remarkable, however, 
that the existing and planned more accurate measure- 
ments can solve these problems and help us to learn how 
the Universe behaved when it was a tiny fraction of a 
second old. 
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